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Resumen 


El presente trabajo proporciona dos gráficas con base en el diseño 
tradicional de una laguna facultativa. En la primera se establece un 
factor para el cálculo del tiempo de retención; en la segunda se 
determina la demanda bioquímica de oxígeno en el efluente del sistema, 
lo anterior con el fin de simplificar el diseño. Se aplicó el criterio del 
método gráfico, luego los valores obtenidos fueron verificados con los de 
la metodología tradicional. Se comprobó la simplificación en el 
dimensionamiento de una laguna facultativa. Es importante indicar que 
el presente estudio considera únicamente el análisis de la materia 


orgánica para el dimensionamiento. 


Palabras clave: laguna facultativa, tiempo de retención, demanda 


bioquímica de oxígeno, diseño gráfico. 


Abstract 


The present paper provides two graphs based in the traditional design of 
a facultative pond. The first one establish a factor for the calculation of 
the retention time; the second one determines the biochemical oxygen 
demand in the effluent of the pond, the aforementioned in order to 
simplify the design. The graphic method was applied, later the data 
obtained were compared with the traditional method. It was 
demonstrated the simplification in the design of a facultative pond. It is 
important to mention that the present paper only considers the organic 
matter for the design. 
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Introducción 


Los costos de mantenimiento y operación de las lagunas de 
estabilización son bajos, por esta razón estos sistemas de tratamiento 
son los más empleados; además, cumplen con la calidad del agua 
residual tratada para reutilizarla en la agricultura (Bracho, Aldana, 
García, € Herrera, 1997; Mara, 2006; Gemitzi, Tsihrintzis, Christou, € 
Petalas, 2007). Una desventaja importante en estos sistemas de 
tratamiento es que se necesitan grandes extensiones de terreno (Wood, 
Greenfield, Howes, Johns, € Keller, 1995; Agunwamba, 2001; Mara, 
2004; Kaya, Dilek, 8 Gókcay, 2007; Naddafi et a/., 2009). Los sistemas 
lagunares tienen como propósito disminuir la concentración de la 
materia orgánica, también conocida con la demanda bioquímica de 
oxígeno (DBO), en donde los procesos pueden ser tanto aerobios como 
anaerobios (Metcalf 8 Eddy Inc., 1991). Según Rolim (2000), las 


lagunas de estabilización se clasifican en anaerobias, facultativas y de 
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maduración, lo anterior en función del contenido de oxígeno. Las 
primeras tienen como característica principal la remoción de grandes 
cantidades de materia orgánica. De acuerdo con la Comisión Nacional 
del Agua (Conagua, 2007a) los estanques anaerobios son profundos (de 
2.5 a 5.0 metros) y la eliminación de la DBO se lleva a cabo por medio 


de bacterias metanogénicas. 


Las lagunas facultativas apoyan su funcionamiento por medio de 
bacterias, algas y la fotosíntesis, en donde la DBO se convierte en algas 
y productos inorgánicos. Una característica importante es que en la 
parte inferior prevalecen condiciones anaerobias; mientras que en la 
parte superior son aerobias (Rolim, 2000). Por otro lado, el viento 
favorece la mezcla, un reparto uniforme de la DBO, bacterias, algas y 
oxígeno disuelto. La consecuencia de lo indicado es un eficaz 
tratamiento (Conagua, 2007b). Sin embargo, estudios experimentales 
han demostrado que tiene más influencia el tamaño de las estructuras 
de entrada que la acción del viento, pues éste es muy variable en su 


dirección y velocidad (Shilton 8 Harrison, 2003). 


La función principal de los estanques de maduración es la 
eliminación de los coliformes fecales. El proceso que se establece es 
aerobio en toda la laguna y el oxígeno es generado por algas y la 
aireación superficial. Este tipo de lagunas son de poca profundidad (de 
0.9 a 1.5 metros) (Metcalf € Eddy Inc., 1991; Conagua, 2007a). De 
acuerdo con Heredia (1985) y Conagua (2007a), los coliformes fecales 
generan enfermedades, como cólera y tifoidea, entre otras. Los métodos 
de diseño que se aplican actualmente están basados en la carga máxima 
por unidad de área (Olukanni, 2011). 
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De los diseños de sistemas lagunares con gráficas 


Swamee y Ojha (1991) reportaron una curva que representa el 
fenómeno de esfuerzo de la DBO como resultado de reacciones 
bioquímicas complejas entre el conjunto heterogéneo de biota y el 
sustrato orgánico complejo. Esta curva muestra una fase de retraso, 
seguida de una curva ascendente terminada en una meseta, y de nuevo 
una curva creciente de tipo saturante, que produce un efecto 
carbonatado de la DBO. Para el sustrato rico en compuestos 
nitrogenados, hay un aumento adicional en la curva de DBO, que satura 
en el valor de la DBO final de la segunda etapa. Ellis y Rodríguez (1995) 
determinaron ecuaciones de regresión múltiple para lagunas facultativas 
y de maduración, donde se relaciona la eliminación de la DBO y 
organismos coliformes fecales a una serie de parámetros ambientales y 
de otros tipos. Kehl, Wichern, Lúbken y Horn (2009) presentaron un 
análisis estadístico para lagunas facultativas aireadas, facultativas y 
lagunas anaerobias por medio de método de análisis Monte Carlo. 
Presentan además algunos análisis estadísticos de los datos, como 
histogramas, gráficos de dispersión y tablas con los valores calculados 
para las variables de las lagunas de estabilización, a partir de los datos 
aleatorias de entrada. Morales, Domínguez y Mata (2016) publicaron un 


estudio estocástico: los autores introducen el factor dentro del modelo 
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de lagunas facultativas por medio del método gráfico. Se aplica la 
gráfica con las isolíneas de fiabilidad obtenidas para el diseño 
estocástico del tratamiento mediante estanques con lemnáceas de los 
efluentes de las lagunas facultativas primarias de los municipios de 


Ranchuelo y Caibarién, en la provincia de Villa Clara, Cuba. 


Del diseño de lagunas facultativas con metodologías 


reducidas 


En la revisión bibliográfica realizada sólo se encontraron dos diseños de 
sistemas lagunares reducidos: a) manual simplificado para el diseño, 
operación y evaluación de lagunas de estabilización de aguas residuales 
de Collí, Rico, Rivas, Escalante y Luyendijk (1992), que considera, en 
términos generales, la calidad del agua residual en la entrada del 
sistema de lagunas, diseño, construcción, operación y mantenimiento; 
b) Treviño-Cansino y Cortés-Martínez (2016) reportaron otro criterio 
simplificado: método de diseño reducido para lagunas de estabilización, 
que estima como valores constantes la evaporación diaria (5 mm/día) y 
la profundidad (1.5 m). Además, es posible definir las dimensiones 
promedio de los taludes. Cabe mencionar que los dos criterios de diseño 


se basan en la metodología tradicional adoptada por la Conagua 
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(2007a). Se presentan las ecuaciones definidas para el último diseño 


exclusivamente de la laguna facultativa. 


Volumen: 
V = BQ; 
Donde: 


E 
Área superficial: 

Asuperior = Bsuperior > Lsuperior 
Ancho superior: 

Bsuperior = € Q;¡ + 1.5 * pendiente 
Longitud superior: 


Lsuperior = 32 Q¡ + 1.5 * pendiente 


(1) 


(2) 


(3) 


(4) 


(5) 
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Donde e es un factor para el cálculo de las dimensiones de la 


laguna que depende de la temperatura. 


Gasto en el efluente: 


Q. = wQ; — 0.03€ Q;¡ * pendiente — 0.01125 * pendiente”2 (6) 


Donde w es un factor para simplificar el cálculo del gasto en el 


efluente. 


Tiempo de retención hidráulico: 


Os = f (7) 


Número de coliformes fecales: 


(8) 


Donde y es un factor que depende de la relación largo-ancho, 


temperatura, y de la constante adimensional (a) y del factor de 
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dispersión (d). Lo anterior para realizar el cálculo de los coliformes en el 


efluente de la laguna. 


DBO en el efluente corregido por evaporación: 


DBO, =TDBO, (2) (9) 


e 


Donde r es un factor para el cálculo de la DBO en el efluente que 
depende del tiempo de retención (0) y de la constante de decaimiento. 
El objetivo del presente trabajo es determinar dos gráficas con base en 
el diseño tradicional de una laguna facultativa. La primera tiene como 
propósito establecer un factor para el cálculo del tiempo de retención; 
mientras que la segunda determina la demanda bioquímica de oxígeno 
en el efluente de un sistema de tratamiento con lagunas de 
estabilización. La aportación del presente estudio es la deducción de un 
método gráfico para el dimensionamiento de una laguna facultativa, a 
partir del método reducido propuesto por Treviño-Cansiño y Cortés- 
Martínez (2016). Una ventaja importante con el criterio presentado es la 
simplificación, rapidez y facilidad para calcular el dimensionamiento de 


una laguna facultativa. 


Materiales y métodos 
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El texto consultado para el diseño fue “Paquetes tecnológicos para el 
tratamiento de excretas y aguas residuales en comunidades rurales” 
(Conagua, 2007a). 


Nomenclatura para el diseño de la laguna facultativa 


C.O. = carga orgánica. 
Q; = caudal en el influente de la laguna (m*/día). 


DBO, = concentración de la demanda bioquímica de oxígeno en el 


influente (mg/l). 

1 000 = factor de conversión. 

A» = Carga orgánica superficial (gDBO5 m*/día). 

T = temperatura mínima media mensual del aire (* C). 

L, = concentración de DBOs en el influente de la laguna (mg/l). 
Af= área en m'. 

Qmea = caudal en el influente en m*/día. 

Z = profundidad en m. 


V = volumen de la laguna (m?). 
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O; = tiempo medio de retención hidráulico, laguna facultativa (días). 
X = relación entre longitud y anchura. 
Bprom = ancho promedio de la laguna en m. 


Lprom = longitud promedio de la laguna en m. 


Bsup = ancho superior en m. 


Lsup = largo superior en m. 


Asup = área superficial en m?. 


Qe = caudal en el efluente de la laguna (m*/día). 


e = evaporación (mm/día). 
d = factor de dispersión adimensional. 
a = constante adimensional. 


N. = coliformes fecales corregidos por evaporación en el efluente de la 
laguna (NMP 100/ml). 


Ni = coliformes fecales en el influente de la laguna (NMP 100/ml). 
Ns/No = número de coliformes fecales en el efluente (NMP 100/ml). 
e = 2.7182818. 

Kf = constante de decaimiento a una temperatura en cualquier día. 
DBO2 = concentración de la DBO); en el efluente de la laguna (mg/l). 
NMamprF = número de mamparas en la laguna facultativa. 


Talud = relación de bordos 2:1. 
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Diseño de laguna facultativa (flujo disperso. Método de 


Yáñez) 
a) Carga orgánica: 
C.0.= Q;¡(DBO;)/1 000 (10) 
b) Carga superficial de diseño: 
As = 250(1.085)7720 (11) 
c) Área de la laguna facultativa: 
As = 10 Li¡Qmea/As (12) 
d) Volumen de la laguna: 
V = (A4)Z) (13) 


e) Tiempo medio de retención hidráulico: 
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Os = V/Q; (14) 


f) Dimensionamiento. Relación largo ancho X = 3: 


Bprom = 2 (15) 


Lprom = al (16) 


g) Para el ancho y largo: 


Bsup = Bprom + (Z)(Talud) (17) 


Lsup = Lprom + (Z)(Talud) (18) 


h) Cálculo del área superficial: 


Asup = (Bsup)(Lsup) (19) 


¡) Gasto en el efluente: 
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Q. = Q; — 0.001 Asup 


j) Remoción de coliformes fecales: 


X = 3 relación largo ancho 


Xx 
— —0.26118+0.25392(X)+1.0136(X)? 


k) Coeficiente de reducción bacteriana: 


Kb = 0.841 (1.075)7720 


|) Constante “a”. Se determina a través de la fórmula: 


a=4 1+ 4K,0fd 


m) Coliformes fecales en el efluente de la laguna facultativa: 


Nf  4aeli-a/2d) 


No (1+a)? i 


n) Coliformes fecales corregidos por evaporación: 


(20) 


(21) 


(22) 


(23) 


(24) 


(25) 
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N¿ = (No)(Q:)/Qe 


ñ) Concentración de la DBO en el efluente de la laguna: 


eN Kf35 
Kf = (1.085)35-T 


(25a) 


(26) 


0) Concentración de la demanda bioquímica de oxígeno en el efluente de 


la laguna: 


DBO; 
KfOf+1 


DBO6+s = 


p) Eficiencia de remoción de la DBO: 


_ (DBO¡- DBO2) 100 
DBO;¡ 


% 


q) DBO corregida por evaporación: 


DBO, = (DBO¡MQ1)/Q6 


(27) 


(28) 


(29) 
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Análisis matemático para el método gráfico de la 


materia orgánica y el tiempo de retención 


Dimensionamiento: tiempo de retención (0). 


Para determinar la gráfica del tiempo de retención se pretende 
realizar una función que tiene de variable independiente la temperatura. 
Para realizar lo anterior, se toma la Ecuación (2), que representa £ del 
método reducido publicado por Treviño-Cansiño y Cortés-Martínez 
(2016). Si se toma la Ecuación (7), que define el tiempo de retención en 
el método de diseño reducido, y se sustituye en la Ecuación (2), se 


tiene: 


z*DBO¡10 


Of = de 


30) 


Por último se sustituye la fórmula de As (Ecuación (11)) en la 


Ecuación (30) y se tiene: 


D0= z*DBO¡10 (31) 


f — 250(1.085)T-20 


Según el método reducido, se fija como constante la profundidad 
(Z) a 1.50 m y a diferencia de éste, se puede calcular el tiempo de 


retención para las diferentes concentraciones de DBO, indicadas por la 
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literatura (Metcalf 8 Eddy Inc., 1991). Se considera la Expresión (26) 
del método de diseño reducido como constante de decaimiento a 
cualquier temperatura, y tomando la Ecuación (31) y expresándola en 


función del tiempo de retención y la DBO; se obtiene: 
10) _ 1.5(10) 
( /DBo,) — 250(1.085)7-20 (32) 


Luego, simplificando la Ecuación (32): 


10) 0.30672276 
DBO;  (1.085)7 (33) 


Graficando con Excel la Ecuación (33) se tiene la Figura 1: 


O/DBOI 


Nino 
/DBO 


/ 


nm 14. >” 35 27 35 20 


TEMPERATURA EN *C 


174 
Tecnología y ciencias del agua, ISSN 2007-2422, 11(2), 158-189. DOI: 10.24850/j-tyca-2020-02-04 


2020, Instituto Mexicano de Tecnología del Agua 
Tecnología y 


CienciaszAgua Open Access bajo la licencia CC BY-NC-SA 4.0 


(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


Figura 1. Gráfica para el cálculo del tiempo de retención. 


Para la DBO en el efluente (DBO,), se toma como base la Fórmula 
(9) de la metodología y se elimina el factor de corrección por 


evaporación (Q;/Q.), y resulta la Expresión (34): 


DBO, =7x* DBO; (34) 


Donde según Treviño-Cansiño y Cortés-Martínez (2016): 


dl zz (35) 


Se realiza la sustitución de la Ecuación (35) en la Ecuación (34) y 


se obtiene la Función (36): 


) x* DBO; (36) 


En la Fórmula (36) se sustituye la constante de decaimiento Kf 
por la Expresión (26) de la metodología tradicional y el tiempo de 


retención O por la Fórmula (31) y se tiene: 


1 
[ Kf35 ll 1.5(10)DBO¡ )+] 
(1.085)35-7 250(1.085)1—20 
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Posteriormente, se simplifica la Fórmula (37) y se sustituye la 


constante adimensional kf35 = 1.2, y resulta la Expresión (38): 


DBO; 
DBO0. = 0.0211780721(DBO¡)+1 (38) 


La Figura 2 muestra la gráfica para determinar la concentración de 


materia orgánica en la salida del sistema de tratamiento. 


DBO! 


Figura 2. Gráfica para calcular la DBO en el efluente de la laguna 


facultativa. 
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Aplicación del método gráfico 


Para la calcular el tiempo de retención, se obtiene el factor (Of/DBO0i) de 


la Figura 2 y se multiplica por la DBOi del caso en estudio: 


0) 


A partir del tiempo de retención se puede deducir el volumen de la 


laguna al multiplicarlo por el gasto en el influente (Qi): 
V =07+*Q; (40) 


El área de la laguna se obtiene al despejar la Fórmula (13) y 
fijando la profundidad Z = 1.50 m: 


ma (41) 


La demanda bioquímica de oxígeno en el efluente (DBO,) se 
obtiene de la Figura 2, con la DBO, como dato conocido. Para obtener el 


gasto en el efluente y las dimensiones de la laguna se aplican las 
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fórmulas del método tradicional (Ecuación (15), Ecuación (16), Ecuación 
(17), Ecuación (18) y Ecuación (20)). 


Ejemplo de aplicación 


Se solicita diseñar una laguna facultativa para una comunidad rural, 
utilizar la metodología tradicional y generar las gráficas para el tiempo 
de retención y la materia orgánica considerando los siguientes datos: 
Gasto Q = 231 m*/día; demanda bioquímica de oxígeno en la entrada del 
sistema de tratamiento DBO, = 220 mg/l. Temperatura T = 11.80*C 
(Conagua, 2007b). Coliformes fecales en la entrada de la laguna: 10 
000 000 NMP/100 ml (Metcalf 8 Eddy Inc., 1991). Finalmente, la 
evaporación con 5 mm/día. 


Resultados 


En la Tabla 1 se muestran los resultados del diseño. 


Tabla 1. Resultados del diseño de una laguna facultativa con el criterio 


tradicional. 
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Resultados 
Datos 
Laguna facultativa 
Q; Nf/No;, | DBO; T Z X As Af Vol. O; 
231 10 000 220 11.8 |1.5 13.00 | 128.06 |3 968.43 | 5952.64 | 25.77 
000 
Resultados, laguna facultativa (continuación) 
d kb a Big Lap Q. Nf/No DBO. Área 
0.31 | 0.46 3.99 39.377 |112.11 | 208.93 | 58 391.01 | 42.98 4 413.87 


Diseño gráfico 


Se utiliza la Figura 1 y Figura 2 para obtener el factor (Of/DBO;) y la 
DBO.,, respectivamente. La Figura 3 y Figura 4 ilustran la aplicación al 
indicar la intersección con la curva. Luego se sigue el procedimiento 
indicado en el apartado de Aplicación del método gráfico: Ecuación (39), 
Ecuación (40) y Ecuación (41). 
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O/DBOI 


0.15 Si ——- 0/DBOi 


0 5 10 15 20 25 30 35 40 
TEMPERATURA EN *C 


Figura 3. Ejemplo para el cálculo del tiempo de retención de una laguna 


facultativa. 


DBOE 


——-DBOe 


10 60 110 160 210 260 310 360 410 460 
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Figura 4. Ejemplo para calcular la DBO en el efluente de una laguna 


facultativa. 


En la Tabla 2 se muestran los resultados del diseño gráfico. 


Tabla 2. Resultados del diseño de laguna facultativa con el método 


gráfico 
Datos Resultados 
Laguna facultativa 

Q; Nf/No; | DBO; T Z X DBO, |0/DBO; Os 

231 10 000 000 | 220 11.8 1.5 3.00 38.87 0.117 25.77 
Resultados, laguna facultativa 

d kb a Bin Lap Q. Nf/No DBO, Área 

0.31 [0.46 3.99 39.37 | 112.11 | 208.93 5 8391.02 | 42.98 4 413.87 


Como se observa en la Tabla 1 y Tabla 2, los resultados son 


similares: el mismo tiempo de retención de 25.77 días; demanda 


bioquímica de oxígeno corregida por evaporación de 42.98 mg/l; área, 4 


413.87 m?; de igual forma las dimensiones. Cabe señalar que el método 


gráfico es más rápido y sencillo de calcular que el tradicional, por lo que 


puede utilizarse de forma confiable. Otra aplicación es para verificar los 


resultados obtenidos con la metodología tradicional; al respecto, Oliveira 


y Sperling (2011) reportaron que cuando los sistemas lagunares no 


cumplen con la calidad del agua residual tratada se debe a errores en su 
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diseño, además de la operación. Ellis y Rodríguez (1995), Kehl et al. 
(2009), Swamee y Ojha (1991) presentaron distintos criterios de diseño 
con diferentes tipos de gráficas; mientras que en el presente estudio 
sólo se consideró el método tradicional de diseño adoptado por la 
Comisión Nacional del Agua y el Instituto Mexicano de Tecnología del 
Agua, pero disminuyendo las expresiones matemáticas a su mínima 
expresión posible, para luego elaborar las gráficas para el tiempo de 
retención y la remoción de la materia orgánica. Una ventaja del estudio 
propuesto es que de las 18 ecuaciones establecidas en el diseño 
tradicional de una laguna facultativa, se redujo a nueve. Las ecuaciones 
suprimidas por el método gráfico fueron las siguientes: (10), (11), (12), 
(13), (14), (26), (27), (28) y (29). Mientras que en el método de diseño 
reducido de Collí et al. (1992) reportaron 10 ecuaciones. La calidad del 
agua residual en el efluente de la laguna facultativa es la NOM-001- 
SEMARNAT-1996, que indica los límites máximos de contaminantes que 
pueden descargarse en los cuerpos receptores; dicha norma fue 
publicada por el Diario Oficial de la Federación (DOF, 1997). Los límites 
máximos permisibles de materia orgánica para uso en riego agrícola que 
indica la citada norma es de 74 mg/l como promedio mensual; por lo 
tanto, la laguna facultativa cumple con la regulación, pues la 


concentración en la salida del estanque resultó de 42.98 mg/l. 


De los coliformes fecales 
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La normatividad señala como límite máximo permisible 1 000 NMP/100 
ml como promedio diario. Los valores de coliformes fecales en el 
efluente del sistema resultaron por arriba del valor permitido, por lo que 
es recomendable el diseño de otra laguna a facultativa o de maduración 
(Rolim, 2000; Conagua, 2007b). Otra alternativa es el uso de mamparas 
o deflectores. Killani y Ogunrombi (1984); Pedahzur, Nasser, Dor, Fattal 
y Shuval (1993); Muttamara y Puetpaiboon (1997); Zanotelli, Medri, 
Belli-Filho, Perdomo y Costa (2002); Shilton y Harrison (2003); 
Sperling, Chernicharo, Soares y Zerbini (2003); Shilton y Mara (2005); 
Abbas, Nasr y Seif (2006); Cortés-Martínez, Treviño-Cansino, Luévanos- 
Rojas, Luévanos-Rojas y Uranga (2014a); Martínez, Cansino, García, 
Kalashnikov y Rojas (2014) reportaron que se incrementa la eliminación 
de contaminantes con el uso de mamparas. En la Figura 5 se muestran 
las dimensiones de la laguna facultativa calculadas con la metodología 


tradicional y con el diseño gráfico. 


Laguna facultativa 
112.11 


39.37 


NF/No; = 106 666 600 NF/No = 58391.02 
DBO; = 220 DBOe = 42.98 
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Figura 5. Ilustración de las dimensiones de la laguna facultativa, 


aplicando ambas metodologías. 


Conclusiones 


De acuerdo con los objetivos planteados, se determinaron dos gráficas 


para el diseño de una laguna facultativa. 


La simplificación en el diseño de la laguna facultativa es de 
consideración, ya que las ecuaciones disminuyeron el 50% 


aproximadamente. 


Con el propósito de simplificar todavía más el diseño se 
recomienda llevar a cabo otro estudio, como trabajo a futuro, donde se 


consideren los coliformes fecales. 


Es importante aclarar que el presente estudio refiere a un análisis 
matemático, por lo que la calidad del agua puede diferir del esperado, 
aunque también influyen la ubicación de entrada y salida de las 
estructuras, la geometría de las lagunas, la operación y. el 


mantenimiento. 
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